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A semiclassical theory based on the Boltzmann transport equation for a two-dimensional electron 
gas modulated along one direction with weak electrostatic or magnetic modulations is proposed. It 
is shown that oscillations of the magnetoresistivity p\ | corresponding to the current driven along the 
modulation lines observed at moderately low magnetic fields, can be explained as classical geometric 
resonances refiecting the commensurability of the period of spatial modulations and the cyclotron 
radius of electrons. 
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Theoretical studies of magnetotransport properties of 
a two-dimensional electron gas (2DEG) spatially mod- 
ulated along a certain direction due to periodic elec- 
trostatic or magnetostatic fields are lasting more than 
a decade. At present the theory of magnetotransport 
in modulated 2DEG is well developed and the most of 
effects observed in such systems at low magnetic fields 
have been explained both by quantum mechanical (in 
a semiclassical limit) P, S S S S IE ^"^^ classical 
[EIEEjOjHI transport calculations, giving consistent 
results. One of a few exceptions is an effect of oscilla- 
tions of the resistivity component pn which corresponds 
to the current driven in parallel with the lines of modu- 
lation. These oscillations were observed along with the 
commensurability oscillations of another resistivity com- 
ponent p±^ corresponding to the current driven across the 
modulation lines, and have the same period and a phase 
opposite to those of pi_. 

The oscillations of p\\ had been immediately explained 
as an effect of a pure quantum nature which originates 
from the quantum oscillations of the electron density of 
states (DOS) in the applied magnetic field B PH- How- 
ever, actual results obtained in these works 0^ Q are 
not consistent with these conclusions. This can be eas- 
ily demonstrated, basing on the results of Ref.0] which 
are reproduced here (Fig. I). It is clearly seen in Fig.l 
that both p± and p\\ exhibit oscillations of the same pe- 
riod at low (unquantizing) magnetic fields where quan- 
tum oscillations of DOS are not resolved, and only effects 
controlled with classical mechanisms could survive. 

It is very well known that periods of quantum oscil- 
lations and semiclassical commensurability oscillations 
(geometric resonance) are different. Their ratio equals 
2t: /\kp where A is the space period of modulations, and 
fci? is the Fermi wave vector of electrons There- 
fore, the observed coincidence on periods of the low-field 
Weiss oscillations of the resistivity component p±_ and the 
weaker anti-phase oscillations of p\ \ gives an evidence that 
for both resistivity components these oscillations are ge- 
ometric oscillations. Such oscillations arise in modulated 
2DEG due to the periodic reproduction of commensura- 
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FIG. 1: Magnetoresistivity components is magnetic field in 
a modulated 2DEG with electrostatic modulations applied 
along the a;-axis 0]. Within the adopted geometry the com- 
ponent pxx, pyy correspond to p±_ and p\\ of the present work, 
respectively. The thik solid line represents the electron DOS 
in units of Do- 



bility between the diameters of electron cyclotron orbits 
and the period of applied modulations A. Consequently, 
we can treat the low-field oscillations of p\ | as an effect ba- 
sically following from classical mechanisms, and develop a 
theory of this effect based on Boltzmann transport equa- 
tion. The results for the magnetic field dependence of p\ \ 
given by the classical transport calculations have to be in 
agreement with the results obtained within the quantum 
mechanical approach in the low-field limit, as well as in 
the case of Weiss oscillations of the magnetoresistivity 
component p^. 

The existing classical theory of magnetotransport in 
modulated 2DEG fails to provide a description of the 
low-field oscillations of p\\, and concludes that this mag- 
netoresistivity component is not affected by the mod- 
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ulations Q|. This gives a motivation for the present 
work. Its purpose is to demonstrate that the most impor- 
tant characteristic features of the low-field oscillations of 
the resistivity pn can be qualitatively reproduced within 
a semiclassical transport theory. We also briefly ana- 
lyze where the earlier semiclassical theory might have 
erred. To simplify the following calculations the effects 
of anisotropy in electron scattering are neglected, and 
the relaxation time approximation is used. It is also 
assumed that the external magnetic field is moderately 
weak so that the electrons cyclotron radius R is consider- 
ably smaller than their mean free path I but larger that 
the period of modulations A, and R >> ^/iX which pro- 
vides preferred conditions for observation of commensu- 
rability oscillations of transport coefficients of the 2DEG. 
To facilitate a comparison of the present results with the 
previous work we intentionally use a theoretical formal- 
ism as close as possible to that adopted in the current 
semiclassical theory 0, 0, . 

We consider first electrostatic modulation with a sin- 
gle harmonic of period A — 271 /g along the y direction 
given by AE{y) = —dV{y)/dy. The screened modulation 
potential V(jj) is parametrized as eV(jj) — eEp sm{gy) 
where Ep is the Fermi energy of the 2DEG. We examine 
weak modulations, so that \egl\ << 1. Electron trans- 
port coefficients for this case were first calculated by 
Beenakker He started from a linearized Boltzmann 
transport equation for the electron distribution function 
^{y, ii) of the form 



+ C[$] = E • V 



(1) 



where ■0 is the angular coordinate of the electron cy- 
clotron orbit; E is a homogeneous electric field applied 
to the 2DEG alongside with the modulating field AE{y); 
the collision term C[^] is written in the relaxation time 
approximation with the relaxation towards the local equi- 
librium distribution, namely: 



1 / 1 /■^'^ 

TV 27r 7o 



(2) 



and the drift term in the case of electrostatic modulations 
applied along the y direction equals 

£»[$] = v{y)smij— + {v' (y) cos i/; + n) — . (3) 

Here, fl is the electron cyclotron frequency and 
the electron velocity vector v(y) has a direction 
u{tlj) = (cos ^, sin?/;), and the magnitude v{y) ~ 
vpy^l + e sin(gy) where vp is the Fermi velocity in the 
unmodulated 2DEG. The Eq.(l) with the collision and 
drift terms of the form (2) and (3) agrees with transport 
equations of later papers [9-11]. The derivation of this 
equation is presented in detail in earlier works (see e.g. 
[10]), and we omit it for brevity. 



The electron current density in the 2DEG modulated 
along the y direction also depends on y and can be writ- 
ten as 

j{y)^Doe^l^ ^v(y,?A)$(y,V'). (4) 

Here, Dq = m/iifi^ is the electron DOS on their Fermi 
surface; m, e are the effective mass and the charge of an 
electron, respectively. 

Following the usual way [10] we write $(y,'0) as: 



$(y,V) = 4>o(V') 



vjy) 

Vp 



- PoT"x(2/, "0) 



(5) 



where po is Drude resistivity, and r is the relaxation time. 
The homogeneous distribution function $o('/') = PoTVojo 
describes the linear response of the 2DEG to the field E 
in the absence of modulations, and the function %(?/, ■(/;) 
satisfies the transport equation: 



D[x] + C[x] = v'{y)v{y)iay. 



(6) 



Here, as well as before, jo is the current density for the 
unmodulated 2DEG. 

We have to remark that these Eqs.(l) and (6) include 
some inaccuracies. Namely, the effect of electrostatic 
modulations is taken into account simply replacing vp 
by viy) which is defined above. This is an oversimpli- 
fication. For systematic consideration we should explic- 
itly introduce into the transport equation an additional 
term which describes the 'electrochemical' field arising 
due to inhomogeneity of electron density caused by mod- 
ulations. As a result the electric field E has to be replaced 
by E — (l/e)V/x where /Lt(j/) is the correction to the chem- 
ical potential of electrons caused by modulations. The 
latter provides the local equilibrium of the 2DEG with 
a stationary electron density {dn/dt = 0). The value of 
this correction fJ-{y) is determined with the continuity 
equation which in this case has the form V • j = 0. 

The above inaccuracies do not change the main approx- 
imation for the distribution function ^(y, tp) but they can 
influence next terms in the expansion of x(y, V") i^i powers 
of a small parameter (fir)"^, as we show below. There- 
fore, using Eq.(6) in further calculations we have to keep 
only those terms in the expansion of xiUj "0) which give 
for the current density results consistent with the conti- 
nuity equation, and all remaining terms in this expansion 
have to be ignored. 

To proceed we expand x(y, ip) in a Fourier series in the 
spatial variable y which gives the system of differential 
equations for the Fourier components. 



dXo 
gRe 



C[xo] 



y-f^t . gRe j uxi 

= — sm^(Xi+X-i)-— cos^l^— , 



dxi , dx-i' 
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= — -^gRmy + — sinV'x±2 

gRt I [dxo , dx±2 
4 \dip dip 



(7) 



The equations determining Fourier components x±n{4') 
for the n > 2 are given by 



± gnR sin 'ijjx±n + 



dx±T 



C[x 



gRe 

cos w 

4 



dtp 
1 — 

dx±n+l , dx±n~l 



gRe . , r , ^ 1 

— Sinip[n{x±n+1 - X±n~l) + X±n+1 + X±r,-lJ 



dip 



dip 



(8) 



Solving these equations, and keeping the terms of the 
order or larger than (egR)^ we arrive at the following 
approximation for the distribution function ip) : 



X(2/, V') = --^^gljlQ {cosip - 
X < sin(gi?cos'0 + gy) — ^^gR 



cos{gR cos ip) cos^ gy — 7; sm{gR cos ip) sin 2gy 



where 



Q 



MgR) 

1- J2(.gi?) 



(9) 



(10) 



and Jo(.gi?) is the Bessel function. Using the obtained 
distribution function we can easily calculate the electron 
current density given by Eq.(4). 

Taking into account only greatest terms in the expan- 
sion of xiUji^) in powers of small parameter (Slr)"^ we 
obtain that only j^, component gets a correction due 
to the modulations along the y direction, whereas the 
component jy remains equal to j° and does not depend 
on spatial coordinates. This agrees with the continuity 
equation V • j = which is necessary to obtain correct 
results for electron transport coefficients in modulated 
2DEG. 

To proceed we need to define effective (averaged over 
the period of modulations) transport coefficients. Here, 
we define the effective conductivity tensor aef / as follows: 



j =<iiy) >= ^ 3iy)dy = (Te//E . (11 



To justify the adopted definition (11) we point out that 
the expressions for transport coefficients obtained either 
with quantum mechanical or classical calculations have to 



be consistent at low magnetic fields. Quantum mechan- 
ical calculations of the magnetoresistivity 0, S IS IE 
give the expression which passes over to a classical con- 
ductivity tensor averaged over the period of modulations. 
Therefore, the latter is an accurate semiclassical analog of 
the conductivity calculated within the proper quantum 
mechanical approach, and our definition of (Jeff agrees 
with that of Ref. [la|. The same definition was used 
previously by Mirlin and Wolfle 0| . 
As a result we obtain: 



^xx _ ^0 , 1 (c-nR\^n- '^O^dR) 



(12) 



where (Tq = 1/po is the Drude conductivity of the elec- 
tron system. Other components of the effective conduc- 
tivity are not influenced due to modulations within our 
approximation. The second term in Eq. (12) arises due 
to modulations. This term exhibits oscillations whose 
origin is explained below. 

It is known that in the presence of both electric and 
magnetic fields centers of electron cyclotron orbits drift 
in parallel with the cross product E x B. The speed of this 
motion Vdr equals E/B which can be easily shown start- 
ing from the equations of motion of an electron subject 
to these fields. This results in one-dimensional diffusion 
of electrons along the direction of the drift. The lat- 
ter brings extra contributions to the current density and 
changes the conductivity of the system. Here, electro- 
static modulation AE(ij) is applied along the y direction 
therefore guiding centers are drifting along the x axis, 
and the electron diffusion caused by the drift reveals it- 
self as a corection to axx component of the magnetocon- 
ductivity tensor. The second term in Eq.(12) represents 
this correction. 

To estimate v^r for the inhomogeneous electric field 
AE{y) we have to average AE{y) over the cyclotron or- 
bit. The mean value is mostly determined with contribu- 
tions from the vicinities of two stationary points at the 
orbit where electron is moving in parallel with the mod- 
ulation equipotential lines (along the x direction), and 
the field AE{y) is varying slowly Correspondingly, 
ti^^ is enhanced when E{y) takes on values of the same 
sign at both stationary points, and reduced otherwise. So 
the diffusion correction in the conductivity cr^ff which is 
proportional to (w^^) ehxibits an oscillatory dependence 
on the magnetic field, an these oscillations are commen- 
surability oscillations due to geometric resonances in the 
guiding center drift. This physical mechanism was pro- 
posed to explain semiclassical Weiss oscillations in the 
magnetoresistivity p± The present result 12] shows 
that the same mechanism controls semiclassical oscilla- 
tions of p|| (see Eqs.(13) ans (14) below). 

The effective magnetoresistivity tensor is defined here 
as the inverse of the effective conductivity introduced by 
Eq.(ll), peff = (^eff- current driven across the 
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modulation lines, the corresponding resistivity is: 

p^=Pyy = Po[l + -^{e9l) 13:^1 ■ (13) 

Assuming that the current flows along the modulation 
lines we get: 



where G{y, ip) satisfies the equation: 



P\\ = Pxx = po<{ 1 - ^{egR)^j 



Mi^X (14) 

For moderately weak magnetic fields {gR >> 1) the ob- 
tained results (13) and (14) describe oscillations of both 
magnetoresistivity components periodic in the inversed 
magnetic field magnitude. The oscillations of p± and p^ 
have the same period in 1/B and the oppositephases 
which corroborates the experimental results of [jj]. The 
amplitude of the oscillations of p|| is considerably smaller 
than that of pj_, and this also agrees with the experiments 
of Ref. 0, and with the results of numerical quantum 
mechanical calculations of Ref. in the limit of a weak 
magnetic field (see Fig.l above). The result for the resis- 
tivity £j_ as we ll ag rees with the corresponding results of 
Refs. |M 1^ Hol lllj obtained within a classical magneto- 
transport theory. 

However, the expression (14) for p|| differs from the 
well known result of the current semiclassical theory. 
To analyze this discrepancy we now calculate the cur- 
rent density taking into account next terms in the ex- 
pansion of the distribution function (9) in powers of 
(fir)"^. Keeping terms no less than (egiJ/Jlr)^ we obtain 
that now all components of the conductivity get grating- 
induced corrections: 



5(7. 



1 + {nr)' 



ra(y); 



-5a,, 



1 



5Uxx 



-5(7rr 



(15) 



Here, a(y) is the y dependent factor of the order of 
{egR)^. It follows from Eq.(15) that both components 
of the current density depend on y and continuity equa- 
tion is violated. On these grounds we conclude that only 
main approximation for xiVj V') can be used in magneto- 
transport calculations, since all remaining terms in the 
expansion of xiUi^) include mistake arising due to orig- 
inal inaccuracies described before. This gives an extra 
justification to our results (13) and (14) for magnetore- 
sistivities. 

At the same time the expressions (15) enable us to 
show the origin of the mistake in the current semiclas- 
sical theory of Refs. ^ ^, 0, ^|. With some formal 
transformations of the transport Eq.(6) we can present 
(q;(?/)) in the form: 

Hy)) = ^l (f(y)sin^G(y,V))# (16) 



D[G] + C[G] 



2v'{y)v{y) 



(17) 



This gives for CTe// and peff expressions absolutel 
consistent with the existing semiclassical theory [T 
Hol Ull . Hence, these expressions could be obtained 
starting from the incorrect expression for the distribution 
function x(j/, ip) which includes terms infiuenced due to 
inaccuracies in the transport Eqs.(l) and (6). This gives 
grounds to seriously doubt the results of earlier works 
[1,0,0,0,01 especially those concerning p\\. 

So, the present analysis shows that simplified transport 
Eq.(l) can be successfully used in calculations of the main 
term in the expansions of transport coefficients in powers 
of (ilr)^^. As for next terms in these expansions, we need 
to modify the transport Eq.(l) (both drift and collision 
terms) to get them. For that we need to carry out con- 
sistent and systematic consideration of effects of internal 
electrochemical field arising due to grating-induced in- 
homogeneity of electron density. This is important for 
redistribution of the electron density at the presence of 
modulations provides the local equilibrium of the system 

Considering magnetic modulations we arrive at similar 
results. Assuming that the magnitude of the magnetic 
field gets a weak spatially dependent correction AB{y) — 
ABsm{gy) {AB/B << 1) we can derive the transport 
equation for the contribution to the distribution function 
of electrons xiUy'^P) which arises due to the modulations: 

D[x] + C[x] = An{y)ivoyf, - voxfy) (18) 

where Af2(y) is the correction to the cyclotron frequency 
due to modulations, and the drift term equals: 

D[x]-^voy^ + {n + An{y))^. (19) 

Writing out the transport Eq.(18) we have taken into ac- 
count that the magnetic modulations do not change the 
magnitude of the electron velocity at the Fermi surface, 
so that v{y) = vp- Due to the same reason the first term 
in the expression (5) now does not depend on spatial co- 
ordinates and equals to the electron distribution function 
for the unmodulated 2DEG. 

Solving Eq.(18) and keeping terms no less than 
{AB / B)'^ in the solution, we calculate the effective con- 
ductivity of the 2DEG at the presence of magnetic modu- 
lations following the way used before to analyze the mag- 
netoresistivities of the 2DEG modulated by a weak elec- 
trostatic field. As a result we obtain that Uxx component 
of the effective conductivity tensor gets a correction due 
to the modulation: 



5<7x, 



( AB 



l-Jl{gK) 



(20) 
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where Ji{gR) is the Bessel function. As in the case 
of electric modulations this correction originates from a 
guiding centers drift. Using Eq.(20) we easily arrive at 
the expressions for the magnetoresistivities p± and p[| for 
the current driven across and along the modulation lines: 

These results describe Weiss commensurability oscilla- 
tions in the limit of weak magnetic fields {gR >> 1). As 
well as for electrostatic modulations, the expression (22) 
reports weaker (due to the small factor [R/lY) oscilla- 
tions of p|| whose phase is opposite with respect to that 
of commensurability oscillations of p^. 

Finally, the novel result of the paper is that it gives 
a qualitative explanation of the low-field oscillations of 
the magnetoresistivity component /9| | in the 2DEG modu- 
lated along one direction within a semiclassical approach. 
It is shown here that these oscillations of p|| at low mag- 
netic fields are commensurability oscillations. At low 
temperatures when the quantum oscillations of the elec- 
tron DOS at the Fermi surface are resolved, Shubnikov-de 
Haas oscillations can be superimposed upon the geomet- 
ric oscillations of the magnetoresistivity. However, this 
does not change a classical nature of the effect itself. 
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